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Fermion creation during preheating in the presence of multiple scalar fields exhibits a range of 
interesting behaviour relevant to estimating post-inflation gravitino abundances. We present non- 
perturbative analysis of fermion production over a 6-dimensional parameter space in an expanding 
background paying particular attention to the interplay between instant preheating (x~^P) an d direct 
fermion preheating (4>-ip). In the broad resonance regime we find that instant fermion production 
is sensitive to suppression of the long wavelength \ modes during inflation. Further, the standard 
scenario of resonant fermionic preheating through inflaton decay can be significantly modified by the 
coupling, and may even lead to a decrease in the number of fermions produced. We explicitly 
include the effects of metric perturbations and demonstrate that they are important at small coupling 
but not at strong coupling, due to the rapid saturation of the Pauli bound. 
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I. INTRODUCTION 
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The evolution of spin 1/2 (and higher) quantum fields 
propagating on curved backgrounds has traditionally 
been considered an esoteric domain of little cosmologi- 
cal interest |lj . In the face of the power and simplicity 
of the inflationary paradigm the added complexities as- 
sociated with their renormalization, regularization and 
solution 1^,0 in the presence of a non-trivial metric g^ 
made them unpopular topics of study. 

The last two years have, however, seen a reversal in 
the fortune of higher-spin fields with the growing belief 
that they may be relevant for cosmology in a number of 
rather profound ways. Decay of the inflaton to massive 
states during inflation Q] can lead to non-trivial features 
in the primordial power spectrum and cosmic microwave 
background (CMB) Q. In inflationary scenarios it is dur- 
ing reheating and preheating - the non-equilibrium end 
to many inflationary models - that fermions first became 
important for the subsequent evolution and nature of the 
universe. 

There are three types of inflaton decay products of 
significant current interest: (i) gravitinos, the supersym- 
metric partner of the graviton present in supergravity 
(SUGRA) ||, (ii) heavy right-handed neutrinos which 
may have mediated lepto-genesis M and (iii) the so-called 



cryptons Q or wimpzillas j^] - massive (~ 10 12 GeV) 
dark matter which may account for the cosmic rays ob- 
served beyond the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off. 

While both heavy right-handed neutrinos and cryp- 
tons provide potential solutions to serious problems, the 
spin-3/2 gravitini are a plague. A typical constraint is 
that their number density 773/2 must satisfy the bound 
n 3/2/ s ^$ 10 where s is the entropy density of the uni- 
verse at reheating for a gravitino of mass 7713/2 ~ 0(100) 
GeV. 

In a perturbative theory of reheating jl(J , this leads to 
the constraint T r < 10 8 GeV [fTlpfl, on the reheat temper- 
ature, T r . Early studies of gravitino production Q fol- 
lowing preheating [|l3-16| argued that thermal gravitino 
abundances were significantly enhanced due to the non- 
perturbative nature of preheating. Nevertheless, it was 
only recently realized that cosmological bounds arising 
from preheating are much tighter than previously imag- 
ined0H2f]g||. 

In these latter works the direct production of gravitini 
was studied through their coupling to the inflaton, qj. 
Nevertheless, in a general SUGRA theory one expects 
many scalar fields <f>i to contribute to the Kahler potential 
K(d>i(j) 1 *). Since the gravitino effective mass is given by 
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where VF is the superpotential and m p i ~ 10 19 GeV is 
the Planck mass, it is natural to ask how other scalar 
fields will affect gravitino production. 

Since we are interested in general results valid also for 
leptogenesis scenarios we will not work directly with the 
gravitino (Rarita-Schwinger) equations but rather with 
the spin 1/2 Dirac equation pij . Fortunately this gives 
useful results for gravitino production due to the na- 
ture of the helicity-3/2 and 1/2 evolution equations. On 
the one hand the helicity-3/2 mode reduces directly to 
a Dirac equation |3|,[l7j and on the other, the danger- 
ous helicity-1/2 mode, which is not suppressed in the 
M p i — > co limit flijf] , behaves like the goldstino in the 
high momentum limit pH-|23|] . 

We will consider a rather general potential leading to 
a fcrmion effective mass 



m e £i = m^ + h X (j) + h 2 x , 



(1.2) 



where is the bare mass of fcrmion, and x are the in- 
flaton and a scalar field coupled to inflaton, respectively. 
The addition of the hi coupling modifies the standard 
production in two interesting cases (i) if (x) ^ during 
inflation or (ii) if x itself undergoes resonant amplifica- 
tion during preheating. We will consider both cases. 

While we will draw general conclusions regarding 
multi-field resonant fermion production t here are certain 
issues we will not address. From Eq. (1.2) it is clear 
that in general the fermion effective mass will have more 
than one natural frequency, analogous to the situation 
with scalar fields. If the natural frequencies (of and x) 
are irrationally related then the effective mass is quasi- 
periodic |p5f . In the scalar case this causes a dramatic 
increase in the strength and breadth of the resonance 
bands. 

Similarly, if the ferm ion is coupled to many scalar fields 
the effective mass (1.2) may become chaotic or stochastic. 
In the scalar case this is known to again enhance the 
resonance if the time correlations are small [|(| . Whether 
these effects occur in the fermion case too is unknown at 
present. 

Instead in this work we will investigate in detail the 
following issues: 

• What is the interplay between direct (through the 
coupling hi) and instant (through hi) p7| fermion pre- 
heating ? 

• How is instant preheating sensitive to the x vacuum 
expectation value (VEV) during inflation ? This will 



strongly alter the fermion effective mass, Eq. (1.2), and 
therefore the nature of stochastic fermion production. 

• How does the inclusion of metric perturbations af- 
fect fermion production ? Are their effects similar to the 
scalar case ? 

To answer these issues requires a 6-dimensional pa- 
rameter space (illustrated in Fig. III). We then study the 



evolution of the spectrum of massive spin- 1/2 fermions 
over this "moduli" space of fermionic preheating. 

In Sec. [0| we discuss this parameter space and present 
the equations for the background geometry and pertur- 
bations. Fermionic preheating is then analysed in the 
massless (Sec. |ll|) and massive cases (Sec. |v|) . 




FIG. 1: Schematic illustration of the potential ( |2.l[ ). 
We consider the inflaton 0, scalar field \ an d fermion ijj. 
Interactions are shown as lines with arrows, leading to a 
6-dimensional parameter space of couplings and masses. 
The new couplings we consider in this paper are shown 
in red and are boxed. 



II. THE MODEL AND BASIC EQUATIONS 



Let us consider the massless chaotic inflation model 
V = ^A0 4 + i 9 W + ffVx 



+ (m + h\(j> + h 2 x)^, 



(2.1) 



where -0 is the fermion field of interest with bare mass 
my, , which is coupled to the inflaton, 0, and a scalar field 
X through the couplings h\ and hi respectively. Typically 
we adopt the value A = 10 -12 for the self-coupling which 
comes from CMB constraints. 

The perturbative (tree-level) decay rates for the infla- 
ton and the x field to fermions are |l(| 
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However, these decay rates cease to be accurate during 
resonance and after a fermion-scalar plasma has formed. 
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Boyanovsky et al. [£8|,g9|, Baacke et al. pO| , pl| and Ram- 
sey et al. |32fl studied this resonant, non-equilibrium pro- 
duction of fermions in the absence of the x held. 

Greene & Kofman |33| gave analytical insight into reso- 
nant fcrmion production and recent work has brought out 
the strong analogy with the scalar stochastic resonance 
regime [|34|-|37| . We extend these works to the multi scalar 
field case where the x ne ld is coupled to both of ip and 4> 
fields. 

We have also included the g 2 <fi 3 x coupling which may 
for example, appear in SUGRA models |58| in addition to 
the usual \g 2 4> 2 x 2 coupling. Due to the existence of this 
term, (x) during inflation and the exponential sup- 
pression for the x ne ld during inflation is avoided |39],fl0| . 
We will see later that g can significantly enhance the 
decays. 

When m.0 = 0, the complete system of equations is 
conformally invariant (modulo the small conformal cou- 
plings needed for <j> and x)- Conformally rescaled fields 
will be capped with a ~ i.e. <p = X = °X? an d 



,3/2 



if) with a the scale factor. 



A. The background and scalar field equations 



In this paper, we work around a flat Fricdmann- 
Lemaitre-Robertson- Walker (FLRW) background with 
perturbations in the longitudinal gauge 



ds z 



-(1 + 2<f>)dt 2 + a 2 (l - 2^)S ij dx i dx j , (2.3) 



where $ and \& are gauge-invariant potentials |4l| , [44| . 
The relation $ = $ follows since the anisotropic stress 
vanishes at linear order fi^j . Decomposing the scalar 
fields into homogeneous parts and fluctuations <p(t, x) — > 
4>(t) + 5<f)(t,x.), x(A x ) — * X(*) + ^x(A x ): one obtains 
the equations for the Fourier modes of the scalar field 
perturbations: 



54> k + ?>H54> k +(^2+ 3A0 2 + .9 V + Qg 2 4>x 
4^6 fc + 2(4, + 3H^ k ~ (2g 2 (b X + 3ffV 2 )5x* 



(2.4) 



S X k + 3HS X k + ( ^2 + 5 2 2 J SXk = 4x$fc 

+ 2(x + 3i7x)*fe - (2g 2 (f>X + Zf^Hk- (2.5) 
A convenient equation for the potential <fr k is 

$*+jre fc =47rG(0ty fc +x*X*), ( 2 -6) 

where G = m~ 2 is Newton's gravitational constant. In 
Eq. ( |2.6| ), we neglect the contribution of fermions since 
it is second order and typically negligible relative to the 
scalar fields due to the Pauli exclusion principle. 



We include the second order energy densities in the 
scalar field fluctuations in the evolution of the back- 
ground quantities via fc-space integrals. The equations 
for the Hubble parameter (H = a /a) and homogeneous 
components of the scalar fields are thus given by Eq] 



H 2 = 



IttG 
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+ -A(0 4 + 60 2 (<5</> 2 ) + 3(50 2 ) 2 ) + -x 2 + -(Sx 2 ) 



+ ^<(V<5x) 2 ) + ^(0V)+3 2 <0 3 x) 



(2.7) 



4> + 3H,p + \4>{<t> 2 +3(6(f> 2 )) 

+ g 2 (x 2 + (Sx 2 m + 3.g 2 (0 2 + W> 2 ))x = 0, (2. 



x + 3H x + g 2 (4> 2 + (H 2 ))x 



(2.9) 



where the expectation values of 5<fi 2 and Sx 2 are defined 



^ J k 2 \5M 2 dk, 



(Sx 2 



2tt 2 



k 2 \S X k\ 2 dk. 



(2.10) 



(2.11) 



As for other definitions of the fluctuational quantities, 
see e.g., Ref. In the present model, both (<5</> 2 ) and 
(Sx 2 ) grow by parametric resonance as shown in Fig. (5), 
and alter the evolution of the background fields. 

We implement backreaction using the Hartree approx- 
imation (i.e. neglecting rescattering effects fl5fl ) and ig- 
nore the backreaction effect due to the production of 
fermions which is typically justified because of the small 
amount of energy in the fermions by the exclusion prin- 
ciple (very massive fermion production being an obvi- 
ous example where this is not necessarily true). The 
Hartree approximation alters the terms 2 and y to 
(£) = (j) 2 + (5<j) 2 ) and (x 2 ) = X 2 + (Sx 2 ) in Eqs. (|J), 
( |2.5| ), respectively. The variances are crucial to ensure 
conservation of energy. 



B. The Dirac equation 

The Dirac equation in curved spacetime is given in 
general by Q 

i-f \d„ + I [7 b , 7 C ] ele cv ,)j J> - m off </> = 0, (2.12) 

where the Latin indices a, b, c, and Greek indices fi, v re- 
fer to the locally inertial (tetrad) coordinates and the 
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general curvilinear (spacetime) coordinates respectively. 
The vierbeins connect these coordinate systems via 



ds 2 = g llv dx tL dx v = r] ab e a e b , 



(2.13) 



where g^ and r] a b are the curved space and Minkowskian 
metrics respectively. The 7 M are the curved-space Dirac 
matrices satisfying the Clifford algebra anti-commutator 
relation = ■ 

Since the spinor ip has no homogeneous component, 
writing the gener al D irac equation to first order in the 
perturbed metric ( |2.3| ) does not yield any terms contain- 
ing <£>fc [in contrast with the scalar case, Eqs. (2.4), (2.5)], 
and simply gives 



(«7 M <9 M - m BS a)ii = , 



(2.14) 



where do denotes the derivative with respect to conformal 
time r\ = J a~ 1 dt. We decompose the ip field into Fourier 
components as 



1 



(2tt) 3 / 2 



d 3k T,s [a s (k)u s (k,T))e 



(2.15) 



ik'Xi 



+ bl(k)v,(Js,f})e 
Imposing the following standard ansatz |33| 
u s {k,r]) = (-i-y^d^ - m eS a)%p k (t)R±(k) 



(2.16) 



where R±{k) are the eigenvectors of the helicity operator, 
which satisfy the relation r y°R±(k) = 1 and k-^ R±(k) — 
±1, we obtain the mode equation for the tp k : 



— ■ + k 2 + (m c ffa) 2 - i— (m e &a) 



d 2 

drj 2 ' drj 
Introducing dimensionless quantities 

A- 2 



4>k = 0. (2.17) 



/ = 



with 



Va^(o) : 



K 2 = 



A</> 2 (0) 



, x = VX(f>(0)r), (2.18) 



m cff = am e s 

= am.^ + hi4> + h 2 x ■ 



(2.19) 



Eq. fl2.17|) becomes 
d 2 



dx 2 ' dx 



which has a WKB-form solution given by 
tpk = a k N + e Jo + /3 k N_e^ J° 



(2.20) 



(2.21) 



where 0| 



2 + f 2 and N ± = l/v/2fi fe (Q fe ±/). The 
comoving number density of produced fermions is given 
in terms of the Bogoliubov coefficients by 




(2.22) 



The initial conditions are chosen to be a k {0) = 1, /3fc(0) = 
0, which corresponds to n k (0) = 0. The Bogoliubov 
coefficients satisfy the relation |afc(i)| 2 + |/3fc(£)| 2 = 1, 
which means that the exclusion principle restricts the 
number density of fermions to below unity, n k {t) < 1. 
In the subsequent sections, we investigate the dynamics 
of fermionic preheating taking into account both the hi 
and hi coupling. 



III. MASSLESS FERMION PRODUCTION 



It was pointed out by several authors [|30 33 that 
massless fermions are resonantly produced by the in- 
teraction between the coherently oscillating inflaton and 
fermion field during preheating. The case of most interest 
qi = h\/\ ^> 1 leads to strongly non-adiabatic evolution 
(f2fc/£! 2 ^> 1) of the ip k only for the short periods when 
m c ff ~ 0. For /12 = 0, this coincides with (j) — and leads 

1/3 

to efficient fermion production for momenta k < q x . 

In most previous works of fermion production, the in- 
teraction between tp and % fields has been neglected. 
However, since the x ne ld can be amplified during pre- 
heating, it is expected that this will lead to the fermion 
production even in the absence of the direct (hi) cou- 
pling. 

The strength of the \ ne ld resonance depends on g. As 
was analytically studied in |4q| , the strongest resonance 
occurs around k — when <p/ A = 2^ 2 where I is an in- 
teger. Hence x particle production is maximally efficient 
even if <7 2 / A = 0(1), which is different from the mas- 
sive inflation model where efficient preheating typically 
requires large couplings g > 3 x 1CP 4 . 



A. light x field: g 2 /X = 0(1) 

Let us first consider fermions with vanishing bare mass 
(to,/, = 0) with g 2 /X = 2 and g = 0, which leads to largest 
growth rate (Floquet index) for the super-Hubble k — > 
modes. We have numerically solved the background and 
perturbed equations and in Fig. the evolution of the 
comoving number density of fermions. Plotted are the 
cases of qi = h 2 /X = 1, 10, 100 and qi = for the mode 
k = 1. We start integrating from the end of inflation, 
and take the initial scalar field values to 0(0) = 0.5m p i 
and x(0) = I0 -9 to p i. Since the x field grows exponen- 
tially through parametric resonance, this changes the fre- 
quency of the fermion field even in the absence of the hi 
coupling. 

In fact, n k increases stochastically until x ~ 50 due to 
the rapid growth of the x field. Not surprisingly we find 
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in Fig. H that increased coupling h 2 leads to larger occu- 
pation numbers. For couplings q 2 > 100, the final num- 
ber density reaches of order unity and further increases 
(<?2 3> 100) cause rik to increase more rapidly. 




FIG. 2: The evolution of the comoving number density of 
fermions versus dimensionless conformal time x during pre- 
heating in the case of 52 = 1, 10, 100 (bottom, middle, top) and 
qi — with g 2 /X — 2 and g = for the mode n = 1. The 
number density grows exponentially until the growth of the \ 
fluctuation stops by backreaction effects. We choose the initial 
values of scalar fields to be 0(0) = 0.5m p i and x(0) = 10 _9 m p i, 
which we use in all figures in which g = 0. Inset: The evolu- 
tion of the inflaton condensate, <j>, vs x. The backreation effect 
becomes relevant for x > 50. 




q 1= 10- 3 , q 2 =0 



For small momentum modes k ~ 0, the excitation of 
fermions is stronger than the sub-Hubble mode in Fig. ^ 
because the nonadiabatic change of the fermion frequency 
is more relevant. After \ particle production stops by 
backreaction effects, the fermion production also termi- 
nates, but does not oscillate. 

We should mention that the final occupation number of 
fermions depends not only on h 2 but also on the value of 
X at the onset of preheating. In the case of g 2 /X = 0(1), 
since the gift term in the lhs of Eq. (2.9) is comparable 
to the Hubble rate, the \ field is hardly damped dur- 
ing inflation j47J. Moreover, even for an initial value 
of x(0) = lO^ripi, which is much smaller than 0(0), 
fermions can be created purely by the h% coupling. 

When the interaction, mediated by h\, between infla- 
ton and fermion is taken into account, fermion produc- 
tion at the initial stage is mainly dominated by the h\ 
coupling unless h 2 3> hi because x is so much smaller 
than (j) at the beginning of preheating. However, para- 
metric amplification of the x field leads to an additional 
growth of the occupation number of fermions. 

In Fig. ||, we plot the evolution of for the cases of 
gi = 10" 5 , q 2 = 10 2 and q x = 10~ 5 , q 2 = with g 2 /\ = 2 
for the mode k = 1. We find that the h 2 coupling assists 
fermion production for x > 30, leading to the final density 
nfc ~ 1. In Fig. ^, the spectrum of the final number 
density is depicted. When q\ — 10 -5 and q 2 — 0, 
does not reach of order unity except for momenta close 
to k = 0. 




q 1= 10- 5 , q 2 =10 2 



q 1= 10- 5 , q 2 =0 



FIG. 3: The evolution of the comoving number density of 
fermions rik during preheating in the cases of qi = KT 5 ,g 2 = 10 2 
(top) and qi = 10 _5 ,g 2 = (bottom) with g 2 /X = 2 and g = 
for the mode « = 1. The initial values of the scalar fields are 
4>(0) = 0.5m p i and x(0) = 10 _9 m p i. The existence of the hi 
coupling clearly assists fermion production in the later phases 
of preheating once the \ field has grown significantly through 
resonance. 



FIG. 4: The spectra of the final comoving number density 
of fermions nt in the cases of (71 = 10 -5 ,q 2 = 10 2 (top) and 
gi = 10~ 5 ,g 2 = (bottom) with g 2 /X = 2 and g = 0. The 
initial values of scalar fields are chosen as 0(0) = 0.5m p i and 
X(0) = 10- 9 m p i. 



In the case of g 2 /X — 0(1), x an d super-Hubble Sxk 
modes are not damped during inflation p7| and have a 
nearly flat spectrum, as opposed to a fc 3 spectrum when 
9 2 /A»l§i 
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This means that when the system enters preheating, 
the growth of \ an d $Xk stimulate the gro wth o f super- 
Hubble modes as can be seen in Eq. ( |2.6D and in 
Fig. ||. The growth rate of the 5\k fluctuation for the 
low momentum mode k ~ is large as /i = 0.2377 for 
.9 2 /A-2@. 

In the standard adiabatic inflation scenario, the 
Bardeen parameter 



H 



xQl) 



(3.1) 



is time-independent in the long- wavelength limit, a re- 
sult which simply follows from energy conservation fl49| . 
Here = 8<j) k + 4>® h /H and Q\ = Sxk + X^k/H are the 
Sasaki-Mukhanov variables. However, in the multi-field 
case where super-Hubble metric perturbations are en- 
hanced, amplification of isocurvature scalar modes leads 
to amplification of on large scales before backreaction 
terminates the resonance as shown in the inset of Fig. |^. 

When 1 < g 2 /X < 3 where super-Hubble modes of the 
$Xk fluctuation are enhanced, super-Hubble metric per- 
turbations are amplified unless x is initially very small. 
In the simulation of Fig. || where the initial value of x 
is chosen as x(0) = 10 _9 m p i, super-Hubble modes 
increase to the order of 0.1. When %(0) > 10~ 8 m p i, we 
numerically find that metric perturbations on large scales 
grow to the nonlinear level $fc ~ 1. 




100 120 



In the rigid spacetime case (found by setting =0), 
the maximum characteristic Floquet exponent for the <5</>fc 
fluctuation is ^ max = 0.03598 in the sub-Hubble range 
3/2 < k 2 /(X(f> 2 ) < -s/3, where <p co is the value when the 
(\> field begins to oscillate coherently. This growth rate is 
one order smaller than that for Sxk- 
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unperturbed metric 



perturbed metric 



K 



FIG. 6: The spectra of the final comoving number density of 
fermions nt in the perturbed and unperturbed ($fc = 0) metric 
with gi=0 and qi = 10 3 for g 2 /X = 2 and g — at x = 50. 
The initial values of scalar fields are chosen as <f>(0) = 0.5m p i 
and x(0) = 10~ 9 m p i. Metric perturbations lead to a decrease 
in the number of fermions produced from \ decays since they 
enhance (5<j> 2 ) which causes the resonance in x to end earlier, at 
smaller values of (x 2 )- 
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FIG. 5: The evolution of the super-Hubble metric perturbation 
$fc with mode k = 10 -26 , and field variances (8\ 2 ) and {8(j> 2 ) 
for g 2 /X = 2 and g = 0. We choose the initial values of scalar 
fields as 0(0) = 0.5m p i and x(0) = 10~ 9 m p i. When the 5\k 
fluctuation is sufficiently amplified, both <I>fc and Sxk begin to 
grow. Finally, backreaction effects of produced particles shut off 
the resonance. Inset: The evolution of the Bardeen parameter 
for the super-Hubble mode k = 10 -26 . The growth of £fc 
is due to the fact that the \ field is light during inflation for 
ff 2 /A = 2@]. 
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FIG. 7: The maximum comoving number density of fermions 
rife with k = 1 achieved during preheating as a function of 
12 = /l|/A. The parameters chosen are g 2 /X — 2,g — 0, 
and qi = with the initial values of <f>(0) = 0.5m p i and 
x(0) = 10 _9 m p i. The fermion occupation number increases al- 
most linearly with q2 until it saturates near nt — 1. The effect 
of the perturbed metric ($£ ^ 0) works to reduce the maximum 
occupation number. 
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In the realistic perturbed metric ^ 0), however, 
the rhs of Eq. ( [j.4| ) leads to the enhancement of the 
S(f>k fluctuation for super-Hubble modes. This makes the 
growth rate of_ (54> 2 ) comparable to that of (8\ 2 ) as is 

is 



X 



If 



10" 



found in Fig. [5]. The final variance 

about two orders of magnitude larger than in the rigid 
spacetime case in the Hartree approximation |o"c| , |5l| ]. 

The rapid increase of the S<pk fluctuation causes back- 
reaction to become dynamically important earlier. This 
means that the final variance of (6x 2 ), which we denote 
((5x 2 )/, is expected to be larger in the rigid case com- 
pared with the perturbed metric case. This is verified 
numerically where we find that (5x 2 ) / is larger in the 
rigid case by about one order of magnitude. 

Clearly this will have an effect on "instant" fermion 
production due to the hi coupling which is modified when 
metric perturbations are taken into account. While the 
spectra of fermions before backreaction sets in are almost 
the same in both the perturbed and unperturbed metric, 
including metric perturbations reduces the maximum oc- 
cupation number of fermions as is evident from Figs. ^|7[ 
although this effect is rather weak. 

This difference practically disappears in the strong 
coupling regime qi 3> 1 due to the saturation of the Pauli 
bound. While we find some difference in the weak cou- 
pling regime, we warn the reader that including second 
order metric perturbations |32| or rescattering effects |53|] 
may change these results, because this is expected to en- 
hance the homogeneous \ field as studied in Ref. fl54[ . 
In addition to this, although we have restricted ourselves 
to the massless inflation, there exist other models of pre- 
heating (3^] where metric perturbations will be amplified. 

As long as the x field is not strongly suppressed during 
inflation and gets amplified during preheating, the basic 
scenario of fermion production will not be changed. Nev- 
ertheless, in order to understand the effect of metric per- 
turbations on fermion production precisely, it is required 
to study the dynamics of fermionic preheating including 
full backreaction effects in broad classes of models, which 
we leave to the future work. 



B. Heavy X field: g 2 /\ > 1 

In the case g 2 /X S> 1, and g = 0, the homogeneous x 
and the super-Hubble 5 X k modes are exponentially sup- 
pressed during inflation since their effective masses are 
much larger than the Hubble rate p4j|35fl . This results 
in the very small values x ^ 10~ 45 m p i at the onset of 
preheating. In this case, fermion production by the ip- 
X interaction is very weak, although the <j)-ip interaction 
leads to standard fermion production. 

In the presence of the g 2 4> 3 X coupling, however, (x) tLP 
|p9| . In f act, neglecting the derivative terms in Eqs. (2.S) 
and (|2.5|), one finds the following relations at the end of 
inflation: 



~ x 2 
9 



- ( - ) 4>, $xk 



-3 |T 2 



(3.2) 



Although these x modes are about (g/g) 2 times smaller 
than those of the inflaton, they undergo parametric am- 
plification during preheating, which leads to fermion pro- 
duction. As an example, we depict in Fig. ||the spectra of 
the final number density of fermions for g 2 /A = 800 and 
g = 10~ 3 g. Since this case corresponds to g 2 /\ — 2l 2 , 
with I — 20, the x modes close to k = are strongly en- 
hanced. In addition to fermion production through the 
hi coupling, the tp-x interaction leads to an extra fermion 
creation as found in Fig. g. 

In the case of g 2 /X ^> 1 with g > 0, the super-Hubble 
5xk fluctuation is about {g/g) 2 times smaller than the 
S(j)k fluctuation, while the sub-Hubble Sxk is the same 
order of 8<j)k- Since the x _ dependent source term on small 
scales in the rhs of Eq. (2.6) is generally larger than the 
large scale modes for g < g, the sub-Hubble metric per- 
turbation grows faster than the the super-Hubble mode 
as shown in Fig. M. 




q 1= 10 z , q 2 =0 



10- 7 



FIG. 8: The final spectra of fermions nt in the cases of 
qi = 10~ 2 ,<72 = 10 2 (top) and 171 = W~ 2 ,q 2 = (bottom) 
with g 2 /\ = 800 and g — I0~ 3 g. The initial values of scalar 
fields are chosen as 0(0) = 0.5m p i and x(0) = -~(g/g) 2 4>(0), val- 
ues whic h we use for all figures in which g 7^ 0, as dictated by 
Eq. ( |3.2| ) We find that the existence of the /12 coupling assists 
fermion production as in the case where g 2 /\ = O(l). 



For g < 10 _2 <7 amplification of the Sxk fluctuation ter- 
minates by backreaction effects before the super-Hubble 
metric perturbation begins to grow. Larger values of g 
would enhance super-Hubble metric perturbations but 
by requiring the inflaton to have a positive effective 
mass, one obtains the following relation making use of 
Eqs. (O) and (O) : 



< 



A 

2^ 



1/4 



(3.3) 
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This means that g/g cannot safely exceed unityQ for 
,9 2 /A 3> 1, which restricts the strong enhancement of 
super-Hubble metric perturbations, and protects the 
CMB from preheating. However, as long as g/g is not 
much smaller than unity and the x held escapes from the 
inflationary suppression, we can say that fermion produc- 
tion is possible only by the /12 coupling as in the case of 
g 2 /\ = (D(l). 
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FIG. 9: The evolution of metric perturbations <&k for the super- 
Hubble mode k = 10 -26 and the sub-Hubble mode k = 1 with 
g 2 /A = 800 and g = 10 _3 (?. The initial values of scalar fields are 
chosen as 0(0) = 0.5m p i and x(0) = — (<?/;?) 2 0(0). In this case, 
the sub-Hubble metric perturbation gets amplified, while the 
super-Hubble mode (and conserved quantity f), remain almost 
constant. 



IV. MASSIVE FERMION PRODUCTION 

A. No % decay: hi ^ 0, hi — 

If the mass of fermion is taken into account, this gener- 
ally works to suppress efficient fermion production. Let 
us first review the case of hi 7^ and /12 = 0. Fermion 
production takes place around the region where the effec- 
tive mass m e ff = m^, + h\(f> vanishes, which corresponds 
to 



This allows an analytical estimate of the maximum 
mass of fermions produced through the hi coupling. Tak- 
ing note that the conformal inflaton field and the scale 
factor evolve as |45l]5l|l 



(4.2) 



= acj> w 4> co cos(O.8472\/A0co7?), 



2ttA 



(4.3) 



3 TOpl 

the effective mass of the fermion field is approximately 



m ef T = m^ + hi 



3 TOpl I 

1 

2ttA (j) co 77 



3(0.8472v / A^co??)- (4.4) 



Since the minimum effective mass is achieved for 
O.8472\/A0co77 = 7T, the upper limit of the fermion mass 
production is 



< TO 

~ 5 



(4.5) 



This indicates that super-heavy fermions whose masses 
are of order 10 17 -10 18 GeV can be produced if we choose 
large couplings with hi close to unity |3(||5| . 

In the vicinity of m c g = where particles are non- 
adiabatically created, we can approximately express m e s 
as 



m cff = hi(f>'^(r/ - 77*), 



(4.6) 



where a prime denotes the derivative with respect to 
rj. Since the expansion rate of the universe is typi- 
cally smaller than the creation rate of fermions around 
TOoff = 0, we can safely set a( f]*) = 1, <x'(?7*) = in 
Eq. fl2.17| ). Substituting Eq. Q) for Eq. 
introducing new variables r 



\fhi¥M - 

y/hx^dri, one obtains (|J3(| 



( pl7| ) and 

- V*), P = 



dr 2 



P 



i>k = 0. 



(4.7) 



The solution of this equation is described by the parabolic 
cylinder function, and the comoving number density of 
fermions is analytically derived as H 



-m^/hi 



(4.1) 



As long as the relation < \hi<j>\ is satisfied, fermions 
are periodically created when the inflaton passes through 
4>*. After the amplitude of inflaton decreases under the 
critical value 0*, fermion particle creation ends. 



*The main concern is that the inflaton must evolve to the 
correct vacuum at <f> = rather then diverging to infinity. 



n k w exp 



7T 2 fc 2 



exp 



\d(f)*/dx\ i/qi 



(4.8) 



This relation means that larger values of qi and \d(j)*/dx\ 
lead to the production of massive particles with higher 
momenta. Massive fermions can be instantly created 
even after one oscillation of inflaton. Their occupation 
number increases to of order unity for low momentum 
modes k ~ 0. Note that if the coupling hi satisfies the 
relation > |fti<ft| at the beginning of preheating, the 



spectrum of Eq. (4.8) can not be applied 
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In order to confirm above analytic estimates, we have 
numerically studied the evolution of the comoving num- 
ber density of fermions. In Fig. [HI we plot the evolution 
of rik for a bare mass ray, = 10 Tr ^ra p i ~ 10 14 GeV with 
g 2 /\ = 2, k = 1, and q x = 10 3 , 10 4 , 10 5 



In this case, the analytic estimate of Eq. (4.5) predicts 
that the coupling hi is required to be hi > 5 x 10~ 5 
for efficient fermion production, which corresponds to 
<7i > 2.5 x 10 3 with A = 1CP 12 . This is elegantly confirmed 
in in Fig. [l^ where fermions are not non-adiabatically 
created for the case of qi = 10 3 . In this case, since the 
bare mass always surpasses the \h\<p\ term during the 
whole stage of preheating, we can not expect stochastic 
increase of n.&. 

On the other hand, for q\ = 10 4 and 10 5 , we find that 
the comoving number density increases to order unity 
when the inflaton passes through 0*. For qi = 10 4 , 
nonadiabatic amplification of fermions occurs only once 
around x ~ 3, after which the amplitude of the infla- 
ton decreases below the critical threshold |</>*|. When 
qi = 10 5 , in contrast, reaches unity several times. If 
one chooses larger values of q\ , nonadiabatic creation of 
fermions takes place periodically for as long as the am- 
plitude of inflaton is larger than |</>*|. 




15 
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FIG. 10: The time evolution of the comoving number density 
of fermions rik with mass my = 10 m p i and the coupling qi — 
10 3 ,10 4 ,10 5 with q 2 = 0, g 2 /\ = 2, and g = for the sub- 
Hubble mode k = 1. The initial values of the scalar fields are 



0(0) = 0.5m p i and x(0) = 10 _9 m 



pi- 



In Fig. we depict the spectrum after one os- 
cillation of the inflaton with bare fermion mass m^, = 
10~ 5 m p i and couplings qi = 10 3 ,10 4 ,10 5 with q 2 = 0, 
g 2 /\ = 2, and g = 0. When qi — 10 4 and 10 5 where 
fermions are created around x ~ 3, the phase spac e den- 
sity of produced particle is described by Eq. (4.8). The 



larger coupling hi leads to the enhancement of higher 
momentum modes, which we can confirm in Fig. [ll]. 

Since nonadiabatic fermion production does not occur 
for qi < 2.5 x 10 3 in this case, the spectrum for qi = 10 3 



does not obeythe analytic estimate of Eq. (fh^) (see the 
inset of Fig. [ll]). Note that the fermion spectra change 
with the passage of time. However, for a short period 
where fermions are created in the vici nity of m e g = 0, 
the spectrum is well described by Eq. (4 



Super-heavy fermions with masses greater than the 
GUT scale ~ 10 16 GeV can be produced if the 
Yukawa-coupling is larger than hi « 5 x 10 -3 by the same 
mechanism described above. These massive fermions 
might have played an important role in leptogenesis sce- 
nario or be candidates of the super-heavy dark matter 
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FIG. 11: The spectra of the comoving number density of 
fermions rik after one oscillation of the inflaton with bare fermion 
mass my = 10 _5 m p i and (?2 = 0. For this mass, simple analyti- 
cal estimates imply a threshold for instant production at around 
qi — 2.5 x 10 3 . The main figure shows nt for qi = 10 4 , 10 s with 
g 2 /X — 2, and g — 0. The spectra match the simple analyti- 
cal result of Eq. (l.i) very well. Inset: vs k for qi below 
the threshold, qi — 10 3 . Essentially no instant production of 
fermions takes place. The initial values of scalar fields are cho- 
sen to be 0(0) = 0.5m p i and x(0) = 10 -9 m p i. 



B. Instant preheating considered: hi = 0, ^0 

Let us next consider the case with hi = and h 2 =/= 0. 
When 1 < g 2 /A < 3, since long-wavelength modes of 
the x field are resonantly amplified, one might assume 
that massive fermions are efficiently produced from the 
beginning of preheating. However, for instant produc- 
tion of s uper-massive fermions, the whole h 2 \ term in 
Eq. (1.2) must be comparable to the bare mass at 
the start of preheating. 

This condition is rather difficult to achieve unless both 
the coupling h 2 and the initial value of \ are large. For 
example, for an initial value x(0) = 10 _9 m p i which we 
used in the simulations of Figs. 0j7l even strong cou- 
plings h 2 ~ 1 yield h 2 x(0) ~ 10^m p i. This indicates 



9 



that super-massive fermions heavier than ~ 10~ 6 7n p i 
cannot be produced at the initial stage. In this case 
however, once the amplitude of the x field grows through 
parametric resonance and the \h%x\ becomes larger than 
m in X ~ production begins to be relevant. 

If g 2 /\ — 0(1) and g — 0, large initial \ values x(0) 3> 
10 m p i will lead to the strong amplification of super- 
Hubble metric perturbations. In order that the model not 
conflict with observations one requires x(0) < 10~ 9 m p i 
which means that instant fermion production of massive 
fermions is rather difficult unless hi is unnaturally large. 

When g 2 /X ^S> 1 with the g coupling, the x field is 
constrained by Eq. ( |3.2[ ). Since the super-Hubble 5\k 
fluctuation is about (g/g) 2 times smaller relative to 5<f>k, 
even the initial value of x(0) ~ 10 _6 m p i does not result 
in the strong enhancement of super-Hubble metric per- 
turbations as shown in Fig. ^. Making use of the relation 
( |3.2[ ) , the effective mass of fermions at the end of inflation 
is approximately given by 



[ ! j) 0(0). 



(4.9) 



When hi and 0(0) are positive, m c s can take negative 
values at the end of inflation for large couplings hi and 
ratio g/g. Since the amplitude of the x field decreases 
at the initial stage of preheating due to cosmic expan- 
sion, instant fermion production takes place when m e g 
vanishes before the x fluctuation is amplified through 
resonance. This case corresponds to 



m p i 



(4.10) 



where we used the value 0(0) = 0.5m p i. 

For example, for a coupling of hi = 1 and g/g = 10 -2 , 
fermions with bare masses of order ra^p ~5x 10 _5 rn p i ~ 
5 x 10 14 GeV are produced at the initial stage. For in- 
stant creation of GUT scale fermions ~ 10 16 GeV, 
we require h 2 (g/g) 2 > 10~ 3 . Due to the existence of the 
g/g term, which is required to be smaller than unity for 
successful inflation, massive fermion production via the 
hi coupling is gen erally hard to realize relative to the hi 
case [c.f. Eq. (|l5l)1. 



It is important to recognise that Eq. (4.10) is the con- 
dition of instant fermion production at the initial stage 
of preheating. Even when hix(0) <C initially, subse- 
quent exponential growth of x leads to fermion produc- 
tion at the later stages of preheating. 

These effects are illustrated in Fig. [12], where we show 
the evolution of n& for fermions with mass m^p — 10~ 6 m p i 
and coupling g 2 /\ = 5.0 x 10 3 , g/g = 3.0 x 10" 3 , ho = 1 
and 0.1 for the mode k = 2. The condition of Eq. ( 4. lCf ) 
is satisfied for hi = 1, which leads to instant fermion 
production found in Fig. [l2[ Since the expansion of the 
universe reduces the amplitude of the x field at the first 
stage of preheating, fermionic preheating is frozen for 



some time once m e ff vanishes at x ~ 2. For low mo- 
mentum modes n ~ 0, increases to unity by instant 
fermion production. 

For x > 30, parametric amplification of the x field 
causes nonadiabatic fermion productio n (se e Fig. fl2"| ). 
For hi — 0.1, where the condition of Eq. ( 4.1CD is not sat- 
isfied, 7ifc does not increase at the initial stage. However, 
rife begins to grow after the x field has been sufficiently 
amplified. Since the homogeneous x field periodically 
changes between positive and negative values, we find 
that fermion productions occurs in the same way even 
for the case of hi < at the final stage of preheating. 
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FIG. 12: Instant production followed by chaotic production 
of fermions. The evolution of the comoving number density of 
fermions nk with bare mass = 10 -6 m p i and the coupling 
h 2 = 1 (solid), hi = 0.1 (dotted) with hi = 0, g 2 /\ = 5.0 x 10 3 , 
and g/g = 3.0 x 10~ 3 for the mode k — 2. The initial conditions 
are 0(0) = 0.5m pl and X (0) = -(g/g) 2 0(O). 
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FIG. 13: The maximum comoving number density of fermions 
nk for the Hubble-scale mode n — 1 as a function of g/g for 
g 2 /\ = 5.0 x 10 3 , hi = 0, and hi = \ with mass = 10 _6 m p i 
and — 10 _8 m p i. As the fermion mass increases, large g/g is 
required before \ — ► ip decays become efficient. 
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In Fig. O, we show the maximum number density 
of fermions as a function of g/g in two cases of = 
10~ 6 to p i and = 10 _8 to p i with g 2 /\ — 5x 10 3 , h\ = 0, 
and hi = 1. We find that large values of g/g are needed as 
one increases the fermion bare mass. In order to produce 
GUT scale fermions m^p ~ 10 16 GeV, we req uire strong 
coupling g/g close to its upper bound of Eq. (3.3). 



H oweve r, the hi coupling makes the second term in 



C. Both couplings taken into account: hi ^ 0, ha ^ 

If both h\ and hi are non-zero, fermionic preheating 
strongly depends both upon the absolute and relative 
strengths of the couplings. In the presence of the g cou- 
pling, the effective mass of fermion at the initial stage of 
preheating is approximately given by 



TOeff « Wty + 



hi - h 2 



~\1 

9_ 

9 



(4.11) 



This relation indicates that the negative coupling hi will 
assist fermion production caused by the positive hi cou- 
pling. Conversely, when hi and h 2 are both positive, 
including the hi coupling generally works as to suppress 
fermion production as long as hi > h 2 (g/g) 2 . 

Let us consider concrete cases for the couplings: hi = 
10" 5 , h 2 = 1, g 2 /X = 5.0 x 10 3 , and g/g = 3.0 x lO" 3 
with niif, — 5.0 x 10~ 7 m p i. When h 2 = 0, fermions are 
stochastically produced at the initial phase of preheating 
when the inflaton passes causes m G ff to vanish. 



1.0 

0.010 
0.00010 

10« 

•A 

10- 8 

10 10 
10 K 
lO 14 

10" 




h 1= 10" 5 , h 2 =1 



15 20 
X 



FIG. 14: /2,2-suppression of instant fermion production. Num- 
ber density of fermions rik vs x in the cases hi = 1CT 5 , h 2 = 1 
(solid), and /ii = 10" 5 , h 2 = (dotted), with g 2 /X = 5.0 x 10 3 , 
g/g = 3.0 x 10 -3 , and = 5.0 x 10~ 7 m p i for the mode 
k = 1. The initial conditions chosen were 0(0) = 0.5m p i and 
x(0) = — (g/g) 2 4>(0). The h 2 coupling reduces the violation of 
the adiabatic condition and removes the initial instant produc- 
tion present when h 2 = for an otherwise equivalent parameter 
set. 



Eq. ( 4.11 ) smaller than the bare mass m^. Hence instant 
fermion production does not take place in this case, al- 
though rtfc increases for x > 30 due to the resonance from 
the h 2 coupling (see Fig. |l4| ). 

This suppression is relevant for large coupling hi rela- 
tive to hi. When hi and hi are of the same order, fermion 
production mainly proceeds by the hi coupling, because 
the ratio g/g is typically smaller than unity. 

When hi = an d h 2 7^ 0, massive fermions which 
satisfy the relation (4.10) are instantly produced. This 
typically requires rather large coupling hi for produc- 
ing massive fermions whose masses are greater than 
rritf, ~ 10 _6 m p i. However, introducing the hi coupling 
even much smaller than hi can potentially alter the dy- 
namics of fermion production. 

For example, in the case of hi = 0, hi = 1, g 2 /X = 
5.0 x 10 3 , and g/g = 3.0 x 10~ 3 with m v , = 10" 5 m p i, 
we can not expect stochastic amplification of fermions at 
the initial stage while increases for x > 40 after the \ 
field is sufficiently produced. Including the hi coupling 
greater than 10~ 4 , fermions are periodically created from 
the beginning of preheating (see Fig. |l5| ). This means 
that the instant preheating scenario discussed in Ref. |27| ] 
can be further strengthened by taking into account the 
hi coupling. 
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FIG. 15: Significant instant production. The time evolution of 
fermion number density for /ii = 10 -3 , h 2 = 1 (solid) and 
hi = 0, hi = 1 (dotted) with g 2 /X = 5.0 x 10 3 , g/g = 3.0 x 10" 3 , 
and — 10 _5 m p i for the mode K = 1. The initial scalar field 
amplitudes are <f>(0) = 0.5m p i and x(0) = ~(3/<?) 2 <^(0). 



V. CONCLUDING REMARKS AND DISCUSSION 

Production of gravitinos and other fermions after infla- 
tion has become an issue of great importance in modern 
cosmology since their creation provides a window into 
the early universe and provides a powerful mechanism 
for distinguishing between different models and regions 
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of parameter space. As such the phenomenon of fermion 
creation deserves to be studied in great detail, especially 
as unusual and non-trivial results are common as the 
complexity of the model scenario is increased. 

With this motivation we have studied spin-1/2 fermion 
production over a 6-dimensional coupling constant pa- 
rameter space with both and without metric perturba- 
tions included. This parameter space spans both the di- 
rect resonant decay of the inflaton into both massive and 
massless fermions and the indirect <f> — > \ i* decays of 
instant preheating. These decays are relevant since they 
model the couplings expected for the gravitino (through 
the Kahler potential) in realistic supergravity theories. 

Previous studies have limited themselves to 3 or 4- 
dimensional parameter spaces in which either there is 
only the cjhip system |3(]|52|,[33| or in which the \ field has 
zero vacuum expectation value during inflation and there 
is no <fr-ip coupling p7j . 

We have verified the general expectation that increas- 
ing the <fi-ip (hi) and x-^> (^2) couplings typically leads to 
enhanced fermion production though with an interesting 
exception. In the case where (x) > and both hi and h 2 
are non-negative, increasing h 2 leads to a decrease in the 
production of fermions due to a weakening of th e deg ree 
of violation of the adiabatic condition, c.f. Eq. (4.11). 

Further, if (x) = during inflation (i.e. g = 0) instant 
preheating is initially sub-dominant relative to direct de- 
cays of the inflaton (Fig. ||) but later becomes important 
due to resonant growth of the x field. If (x) 7^ during 
inflation an initial burst of instant preheating is typically 
followed by a quiescent period and then a sudden chaotic 
burst of fermion production when the x field becomes 
large due to resonance (see Figs. |l^ , |l4| ). 

In addition we found that metric perturbations had an 
important impact on n& at weak coupling (small hi,h 2 ) 
while their effect is weak at strong coupling due to the 
saturation caused by the Pauli exclusion principle. In the 
cases we studied (Figs. |6[ ^) the effect of metric perturba- 
tions was to reduce the production of fermions. This can 
be understood as due to the enhancement of the inflaton 
resonance and hence inflaton variance (5<ft 2 ) when metric 
perturbations are included p4 47 1. This causes backreac- 
tion to end the resonance earlier, before the inflaton and 
X fields have decayed significantly to fermions. 

Fermion production via the h 2 coupling is sensitive 
to the value of x a t the end of inflation. Although we 
considered the massless inflaton model where the x field 
may escape suppression during inflation for g 2 /A = 0(1), 
there exists other scenarios of inflation in which the ex- 
ponential suppression can be avoided. One such scenario 
is a nonminimal coupled x field [^6 47 1. During inflation, 
the homogeneous and super-Hubble x modes exhibit ex- 
ponential growth for negative nonminimal coupling £ as 
pointed out in Ref. |57f| , which may drastically change 
the scenario of instant fermionic preheating described in 
this paper. This amplification of the x field is expected 
to work even for small values of |£| smaller than unity in 
a model-independent way. Nevertheless this requires a 



detailed study including the evolution of metric pertur- 
bations during inflation and preheating, which we leave 
to the future work. 
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